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ASYMPTOTICS AND EXACT FORMULAS FOR ZAGIER 

POLYNOMIALS 

ATUL DIXIT, M. LAWRENCE GLASSER, VICTOR H. MOLL, AND CHRISTOPHE VIGNAT 


Abstract. In 1998 Don Zagier introduced the modified Bernoulli numbers S* 
and showed that they satisfy amusing variants of some properties of Bernoulli 
numbers. In particular, he studied the asymptotic behavior of and also 

obtained an exact formula for them, the motivation for which came from the 
representation of B 2 n in terms of the Riemann zeta function (^{2n). The modi¬ 
fied Bernoulli numbers were recently generalized to Zagier polynomials B* (x). 
For 0 < X < 1, an exact formula for B 2 ^{x) involving infinite series of Bessel 
function of the second kind and Chebyshev polynomials, that yields Zagier’s 
formula in a limiting case, is established here. Such series arise in diffrac¬ 
tion theory. An analogous formula for BJ,^_j_^(x) is also presented. The 6- 
periodicity of BJ^_|_^ is deduced as a limiting case of it. These formulas are 
reminiscent of the Fourier expansions of Bernoulli polynomials. Some new re¬ 
sults, for example, the one yielding the derivative of the Bessel function of the 
first kind with respect to its order as the Fourier coefficient of a function involv¬ 
ing Chebyshev polynomials, are obtained in the course of proving these exact 
formulas. The asymptotic behavior of Zagier polynomials is also derived from 
them. Finally, a Zagier-type exact formula is obtained for BJ^ (“f) + ^ 2 n- 


1. Introduction 


The Bernoulli polynomials are defined by means of the generating function 

3 P. 3] 

Bn{x)z^ 


(|z| < 27r). 


n—0 


In 


( 1 . 1 ) 


Don Zagier introduced the modified Bernoulli numbers 


K :=E 


r—0 


Br 


2r J n + r 


{n G N), 


where Rn(0) are the Bernoulli numbers. In this paper, he proved three 

remarkable results for the sequence 


(A) The value of R* for n odd is 6-periodic: more precisely, it is given by 


n mod 12 
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(B) The fractional part of the number B 2 n ■= — B 2 n satisfies 

B 2 n = - mod 1 (n G N). 

p 

(p + l)|2n 

p prime 


(C) Let n G N. Then is asymptotically equal to (—1)" ^(27r) ^”(2n — 1)! for 
n large, and is given much more precisely by the approximation 

(1.2) Ri (-l)"7ry2n(47r) {n^oo), 

where Yn{z) denotes the Bessel function of the second kind of integer order n defined 
by [27l p. 64] 

(1.3) r„(z) = lim Y^{z), 

v—¥n 


where Y^{z) is the Bessel function of the second kind of non-integer order v defined 

by 

V cos {vtt) - J-v{,z) 

‘ ( \ ’ 

Sin 

with Jjj{z) being the Bessel function of the first kind m p- 40] 


(1.4) 


m=0 


(-l)’”(z/2)2'"+'^ 

m!r(m -I- 1 -I- 1 ^) 


The sign Ri in (HU means that as n —>■ oo, the relative error between the two sides 
decays more rapidly than any polynomial power n~^. In particular, it yields the 
asymptotic 


(1.5) B 2 „--(-l)”7rl2n(47r) (n-)-oo). 

The result on the 6-periodicity of B* for odd n first arose in Zagier’s work p5] . 
where he obtained a new proof, based on the theory of periods of modular forms, 
of the Eichler-Selberg trace formula for the trace of Hecke operator acting on 
modular forms on SL 2 (Z). The method of the proof gave a formula for these traces 
in a form different than the usual, and involved Bernoulli numbers. The special case 
i=\ gave the dimension of Mk (SL 2 (Z)) (fc even) in terms of and the equality 
of this formula with the standard dimension formula required the 6-periodicity to 
hold. 

The result in (^4) was extended in [7] to the so-called Zagier polynomials 



Br{x) 
n + r 


(n > 0), 


where B^{0) = B*. It should be mentioned here that these polynomials were 
briefly studied by Zagier himself, and appear in an exercise in [5l p. 122, Exercise 
20]. Specifically, he obtained some results associated with their generating function 
EZiKi x)z’^, viewed as a formal power series. However, an explicit formula for 
this generating function in terms of the digamma function was first obtained in [7J 
Theorem 3.1]. 

As shown in [3 Theorem 1.2], {i? 2 n-i-i(i)}n(=o is periodic and non-constant pre¬ 
cisely when j G {—3, —2, —1, 0}. The period is 6 for j = —3, 0, and 2 for j = 
—2, —1. Moreover {i3|„_|_^(j)-|-n})]TQ is periodic if and only if j = —4. Equivalently, 
using symmetry result for Zagier polynomials (see (18.311 below), {—B 2 rt+i(j)+^}rj°=o 
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is periodic if and only if j = 1. The only other value of j that yields a periodic exam¬ 
ple is j = —3/2, in which case we get the vanishing sequence since i32n+i(~3/2) = 0. 
Similar periodicity results hold for the sequence — j) — see 

m- 

The second paper in this series [4] was inspired from property (B) above, and the 
arithmetic nature of Zagier polynomials with integer arguments was studied there. 
Let an and 7 „j respectively be the denominators of i?* and B*(j) in reduced forms, 
where / G Z and n G N. Again, two different behaviors, one for n even and another 
for n odd, are observed. The fact that a 2 n+i = 4 for all n G N was proved in [7]. 
It was also conjectured there that 4 divides a 2 n for all n G N. This conjecture was 
established in [3] in the form of the following theorem, which also implies that 
is independent of j. 

Theorem 1.1. Letp be a prime and let Vp{(.) denote thep-adie valuation of I (that 
is, the highest power of p that divides i). For n G N, 

{ 1 ifn = 6 mod 12, 

2 ifn = 0 mod 12, 

0 otherwise. 


The result in (C) motivates the present paper. Zagier 
ogous to the asymptotic formula [U p. 267] 

2(27r)-2-(-l)-+i(2n)! (n 


obtained (IE3, anal- 




2n 


satisfied by the Bernoulli numbers. The above asymptotic formula for Bernoulli 
numbers can be replaced by the following exact formula [a p. 266], m p- 5, 
Equation (1.14)] 


( 1 . 6 ) 


B2n — 


2(-i)»+i(2n)!C(2n) 

(27r)2" 


where ('(s) is the Riemann zeta function. Zagier, in looking for a corresponding 
exact formula for B^n, obtained the following beautiful result [30] 


(1.7) = -„ + ^(^(-l)"7ry2„(4T^) + ^)-icQ) 

, ^ 1 ( yi+i-vi 

^£(£ + 4 ) 2 



The heuristic behind the discovery of this formula is interesting, and the interested 
reader is referred to [30]. The proof of (US is pretty and involves meticulous 
manipulations of the associated infinite series and integrals. 

Note that the companion formula 


( 1 . 8 ) 


B. 


1 


2n+l 


-4 


4 V 2n-f 1 


1 / -3 \ 
”^2 V 2n -h 1 y ’ 


where (^) is the Jacobi symbol, is elementary and has also been established in [30] . 
It can also be rephrased in the form [7] Corollary 10.6] 


(1.9) 


B 


* _ 


(-ir 

4 



^ (2n J- l)7r^ 
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For 0 < x < 1 and n > 1, the Fourier expansion of the even-indexed Bernoulli 
polynomials is given by [22l p. 5] 

( 1 . 10 ) 

m—1 ^ ^ 

which gives (HU as a special case when x = 0. Similarly, for 0 < x < 1 when n > 0, 
and for 0 < x < 1 when n = 0, the Fourier expansion of the odd-indexed Bernoulli 
polynomials is given by |221 p. 5]0 


( 1 . 11 ) 


B2n+i{x) = 2{-ir+\2n+l)lJ2 


sin 27rmx 


(27rTO)^"+i 


The Fourier expansions in (11.101) and (11.111) now raise two natural questions: does 
there exist a generalization of (HU for the Zagier polynomials and, is there 

an analogue of such a generalization for B^^+iix)? 

The primary goal of this paper is to answer these two questions in the affirmative. 
The generalization of HU for the polynomial B^^ix) when 0 < x < 1 is given first. 


Theorem 1.2. Let 0 < x < 1 and n S N. Define 


( 1 . 12 ) 


9{y,r,x) 


jy + l + x- 1 + x){y + 3 +1^ 

Viy- I + x){y + 3 + x) 


Let Yn{z) be defined in HU: and denote by Un{x) the Chebyshev polynomial of the 
second kind. Then, 


BLix) = (-l)"7r^ (F2n(47rm) cos(27rmx) 

m—1 

(^) + (I) + (^) + (^)) 

^ / OO CX3 \ 

22n-\-l g{m, n, x) -I- ^ g{m, n, 1 - x) j . 

\m—l m—1 / 

The analogue of Theorem 11.21 for S 2 n-i-i(^) given next. 

Theorem 1.3. Let 0 < x < 1 and n G N. Let g{y, r, x), T„(z) and Un{x) be defined 
as before. Then, 


OO 

B^n+ii^) = (-l)"7r l2n-n(47rm) sin(27rTOx) 

m—1 



^There is a typo in the version of this formula given in the book. The power of —1 there should 
be n and not n + 1. 
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Remark 1. The fact that the Bessel function series appearing in Theorems 11.21 
and 1 1.31 converge conditionally can be proved as follows. The asymptotic expansion 
of Y^{z) as \z\ ^ oo (see (|2.1I) below) implies that as m —>■ oo, 

(1.13) (-l)"7ry2«(47rTO) -- 

Zy/m 


Now it is well-known m p. 257] that the series ® converges con¬ 
ditionally for Re(s) > 0, which means that each of the series J2m=i and 

sin( 2 ^a:) converges conditionally, and hence the series in Theorems 11.21 
and fOl as well. 

Remark 2. Theorems 11.21 and 11.31 combined with the result [3 Lemma 10.2] 

+ 1 ) = B^{x) + 2^ri-l (2 ’ 

or more generally with 


(1.14) 


i?:(x+fc) = i?:(x)-f 

1=1 



{k e z), 


give exact formulas for B^(x) for any non-integer real values of x. 

In a completely different context of diffraction theory, V. Twersky |261 Equations 
(40), (41)] (see also C. M. Linton [14], [HI Equations (47), (49)]) has obtained one 
of the intermediate results in the proof of our Theorem 11.21 namely (I3.29|) (and 
likewise the corresponding equation occuring in the proof of our Theorem (Oj), 
albeit these are expressed in forms much different than ours, especially since they 
are phrased using the terminology of diffraction grating. However, our proofs of 
these intermediate results are new and completely different from his. Our proofs 
also give new and important results along the way, for e.g.. Lemmas 13.31 and 13.41 

As hinted above, infinite series involving Bessel functions and trigonometric func¬ 
tions are often encountered in studies on the theory of diffraction m, m, [SI, 113 
Appendix D], [TBj. Indeed, the series in the above theorems have arisen [53], [5S] 
in the analysis of the scattering of a plane wave on a diffraction grating with an 
arbitrary angle of incidence. The Schlomilch or Schlomilch type series, such as the 
one in Theorems ll.2l and ll.31 or the ones studied in [2], are conditionally convergent 
and converge extremely slowly. Their partial sums are highly oscillatory. Thus the 
advantage of alternate representations for them, such as the ones in Theorems 11.21 
and fL3l is that they are useful for their fast computation. 

W. v. Ignatowsky m Section 6] studied the series Hm=i Y 2 n{mD), where n S 
N,D > 0 and D is not an integral multiple of 2tt. His case obviously does not cover 
the series in Zagier’s formula ([12]). However, we show that Zagier’s formula can 
be obtained from Theorem O in the limiting case x —>■ 1, the proof of which is 
interesting in itself. Also, it is shown that the limiting case x —>■ 0 of Theorem 11.31 
gives the curious 6-periodicity of i? 2 Tt+i mentioned in property (A), and which is 
equivalent to (11.811 and (11.911 . 

K. Dilcher [6] Corollary 1] showed that the sequence of Bernoulli polynomials 
converges uniformly on compact subsets of C to the sine or cosine functions. This 
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implies, in particular, the asymptotic formulas 

B 2 n{x) ^ 2(27r)“^"(2n)!(—1)"“^ cos(27ra;), 

B 2 n+iix) 2(27r)“(^"'+^)(2n + sin(27ra:), 

for real x as n —>■ oo. 

As an application of our Theorems ll.2l and ll.31 the following asymptotic relations 
for Zagier polynomials are obtained here. 

Corollary 1.4. Let 0 < x < 1. For x ^ as n —>■ oo, 

(1.15) B2n(x) -- (-l)”7rr2„(47r)cos(27rx), 

and 

BL (J) (-l)"+VF2n(8T) ^ i?*„ (J) . 

Also for X ^ as n ^ oo, 

B 2 n+iix) ^ (-l)"7ry2n+i(47r)sin(27rx). 

Lastly, we obtain a Zagier-type formula linking an infinite series involving Y 2 n( 8 TTm) 
with (“i) ^ 2 n- This formula is similar in flavor to Zagier’s formula (11.71) . 

The genesis of this formula is explained in Section [5] 


Theorem 1.5. Let [/„(x) be defined as before. The identity 


-'2n 


-0+«L = 2y:(,-i)vw8™) + ^) 

m—1 ^ V / 


is true for oZZ n G N. 


1 


- n - - U2n-1 7 + U2n-l T ” X 


1 . / I 


72^ 


/ _\ (i 

^ o° ( m -f 4 — ^Jm{m + 8) j 

24"-i ^m{m + 8) 


m—1 


This paper is organized as follows. The preliminary results are collected in 
Section[2j Theorem 1 1.2 1 is proved in Section[3l We do not give the proof of Theorem 
11.31 since the approach is similar to that of Theorem 11.21 Section |4] is devoted to 
deriving Zagier’s formula (11.711 as a special case of Theorem 11.21 Similarly, Section 
[S] contains proof of the 6-periodicity of i? 2 n-i-i resulting from Theorem 11.31 The 
asymptotic properties of Zagier polynomials are proved in Section [51 The Zagier- 
type exact formula for (~|) + ^ 2 n derived in Section [T] Finally the paper 
concludes with three open problems discussed in Section |51 


2. Preliminaries 


The Bessel functions Jv{z), and Yi,{z) have the following asymptotic expansions 
for \z\ —>■ 00 and | a.rgz\ < n [m p. 199]: 


TTZ 


COS It; 


E- 

n—O 


(-ir{n,2n) 


i 2 z) 


2n 


— Sin It; 


E 

n—O 


(-l)"(z/, 2n-f 1) 
(2z)2"+i 
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(2.1) Y,(z) 


1 

2 \ 2 
TTZ 


sin«;Vizip^+cosR;V 2n + 1) 


( 2 zr 


(2z)2"+i 


n=0 ^ ^ n—0 

Here w = z - ^tti^ - and (u, n) = r(nlmJ'-nl\/ 2 ) ■ 

Also, if u —>■ oo through positive real values, then the asymptotic expansions of 
J^(z) and Y^(z) for a non-zero fixed z are given by [20l p. 231] 

1 (ez\^ 


Mz) 


(£)' 


\/ 2 t:v 

2 /ez\-'' 


( 2 . 2 ) YAz)^-J—(^\ 

V TTU V2l// 

The property [501 P- 222, 10.4.1] 

J.n{z) = i-irMz) 

for n G N is used throughout the paper without mention. So will be the facts 

Jo(0) = l, J.(0)=0, 

for Re(u) > 0. 

The Chebyshev polynomials of the first and second kinds are respectively defined 
for n > 0 by the Binet formulas 


(2.3) 


Tnix) = 
Unix) = 


(x + y/x^ — 1)" -I- (a; — y/x^ — 1)” 

2 ’ 

(x + y/x^ — 1)"+^ — {x — y/x^ — 1)"+^ 


2 y/x‘^ - 1 


They are alternatively given by 


r„(cos6>) = cos(n0), UnicosO) = 


sin((n -I- 1)0) 
sin0 


The Poisson summation formula [501 p. 60-61] states that if /(t) is continuous and 
of bounded variation on [0, oo), and if f{t) dt exists, then 

nOO nOO 

(2.4) /(O) + 2 V f{m) = 2 f{t) dt -I- 4 V / f{t) cos(27rmt) dt. 

_ n Jo _ n Jo 


m—1 


m—1'' 


3. Proof of Theorem 11.21 


The proof of (II.7L as given in [30], begins with the representation (11.11) and 
uses the fact that all odd-indexed Bernoulli numbers, except the first one, vanish 
and that the even-indexed ones, namely i? 2 r, can be expressed in terms of the 
Riemann zeta function C(2r). This procedure does not extend very well to the case 
of Bernoulli polynomials, simply because the odd-indexed Bernoulli polynomials do 
not vanish, and leads to two complicated terminating 5 T 2 hypergeometric functions 
which do not seem to produce anything like dnD. 

The idea is to start with the formula 


mni^) = E(-i) 

r—O 


n+r 


I n r 
V 2 r 


B 2 rix) 


+ U2n-1 


(f) 


+ U 2 n -1 


X +1 
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established in [U Theorem 10.1], Observe that only even-indexed Bernoulli poly¬ 
nomials appear in this representation. Separating the term r = 0 and then using 
(ll.lOp yields 

(3.1) = t^+A{n,x) + i (|) + (^)) 

with A{n, x) defined by 


n—1 


A(n,x) = (—^ cos(27rma;) ^ 


(n -I- r)! 


1 


(n — r — 1)! (271771,)^’’+^ 


Upon replacing r by n — 1 — r in the above sum, the function A{n, x) can be 
written in the form 


CXD 


(3.2) A{n,x) = (—^ S'2„(47rTO) cos(27rma;) 

m—1 

where Sn{x) is the Schlafli polynomial [33 P- 285] defined by So{z) = 0 and 


n—a 

2 


O / ^ {n-r-iy. 

= 2^-M- 


r—0 




with 



for n even, 
for n odd. 


The proof is now broken down into a series of lemmas for an easy perusal. A 
new expression for A(n, x) is first presented. 


Lemma 3.1. The function A(n,x) is given by 


OO 

A{n,x) = (— l)"7r ^ U2n(47rTO) cos(27rTOa;) 

m—1 

OO 

+(-!)"« E {2(7 -I- log(27rm)) J2„(47rm) -I- P2n(47rm) — 2 ( 52 n( 47 rTO)} cos(27rma;), 

m—1 


where 7 is Euler’s constant, Yn{x) and Jn(x) are defined bv (II.3p and (11.41) re¬ 
spectively, and the functions P„, Qn are given by [33 p. 341] □ 


and 


Pu{z) 


n—1 






00 

e=o 


z N {fpin -\- i -\-1) — -\-1)) 

2 ) £! {n + e)\ 


Qn{z^ 


{fijn + e + 1) + j) 
he ^ V 2 J A.{n + e)\ 


^We have used Pn{z) and Qniz), instead of the standard notation Tn{z) and Un{z), in order 
to avoid any possible confusion with the notation for the Chebyshev polynomials of the first and 
second kind. 
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where 'ip{z) is the logarithmie derivative of the gamma funetion r( 0 ). 
Proof. For n £ N, the identity [27l p. 340, 10.6] 


n —1 


7ry„(z) = 


(n — r — 1 )! f 


r—0 


(I)' 


+ 


e=o 


I /2 \ (2 log(z/ 2 ) — ip(i + 1 ) — i/;(n + £ + l) 


£\ in + i)\ 

can be written in the form [m p. 340, 10.6(1)] 

(3.3) Sn{z) = -TtY„{z) + 2(7 + log(z/2)) J„(z) + P„(z) - 2Qn{z). 

Now substitute (|3.3I1 . with 2 : = 47 rm and n replaced by 2n, in (13.2|) to produce the 
result. That the right-hand side can be written as the sum of two series follows 
from the first remark after the statement of Theorem 11.31 □ 


The next lemma gives yet another expression for the function A(n,x). 
Lemma 3.2. The function A(ri,x) is given by 

00 

A{n,x) = (-l)"7r E l2n(47rm) cos(27rTOx) 

m—1 

+ (—1)"+^ ^2 ( P 2 n( 47 rm) -I- 2 — Jy( 47 rm) j cos{2'Kmx). 

Proof. The functions Pn{z) and Qn{z) are also given by [13 p. 344] 


~ dn+2k{z'). 


(3.4) 

Pn{z) 

00 ^ 

— ^ ^ ~i^{Jn-\- 2 k{^) Jn- 



k^l 

(3.5) 

Quiz) 

n -| 00 / 

= ^"WE-+E^ 


k(n -I- k) 

Also, for a non-negative integer u, the derivative of Jv{z) with respect to its order 
is given by [13 p. 53, formula (39)] 

(3.6) = (log (I) - + 1)) i{=) - f; + ■'> 


fc=i 


k{k + v) 


-J2k+v{z), 


Now (1531) and (13.61) together give 


1 d 

^ - -b log (271771) - 7/1(277 -b 1)J - —Ji, (471777) 

SO that 

2 (7 -b log (271777)) J2„ (471777) -b P2n(471777) - 2Q2n{^'Xm) 

= 2^7 -b 7/1(277 + 1) - E (471777) -b P2„(47r7?7) -b 2 —J^{4:TTm) 


h'—2n 


The identity 


2 n 


7 / 1 ( 277 -b 1) = -7 + ^- 


i=i 









10 


A. DIXIT, M. L. GLASSER, V. H. MOLL, AND C. VIGNAT 


that appears as entry 8.365.4 in [9] and Lemma [3.II now complete the proof. □ 


The next task is to find an almost closed-form expression for the second series 
in the above lemma. The following two new lemmas, interesting in their own right, 
show that the functions P2ra(47rm) and arise as Fourier coefficients 

in the Fourier expansions of some relatively simple functions. 


Lemma 3.3. For n £ N and 0 < a; < 1, the identity 

OO 

P2n(47rTO) cos(27rTOa;) 

'a: -I- 1 


1 (-ir 

2 n 2?! 


COS ^ C/2n-l (I) +COS ^ 


• COS 


1 — X 


U2n-1 


1 — X 


■ COS 


U2n-1 

2 — X 


X -I- 1 
2 

U2n-1 


2 — X 


holds. 


Proof. Let f{x,n) denote the right-hand side of the above identity. Note that 
/(I — x,n) = f{x, n). The periodization of /(x, n), as a function of x, based on its 
values in [0,1) makes it an even function of x. Hence its Fourier series is given by 

OO 

Oq + Ora COs(27rTOX), 

m—1 


where 


= / f{x,n)dx, Om = 2 / f {x, n) cos{2nmx) dx, 

Jo Jo 

for m > 1. These Fourier coefficients are now computed. Note that the change of 
variable x = 2 cos 0 yields 

t /2 


Vin-X (t 

CN _i/a 

T cos (- 

dx = 2 [ 

/o 

IJ \ i 

J 7T 


Similarly, combining all such integrals, it is seen that 


1 


ao = 


2{-iy 


^7r/2 


2n 
= 0 . 


dsin(2n0) dO 


Also, a similar procedure gives 
(3.7) 

cos(27rmx) 4(—1)" 


Om. — 2 


4(-l)^ 


2n 

fix 12. 


■ dx + 


0sin(2n6l) cos (471771 cos 0) dd 


0sin(2770)| Jo(47r777) + 2^2{—iy J2j{4:Trm) cos{2j6) j d9, 


where in the last step, formula 10.12.3 in EOl p. 226] was used. Now 
l)”Jo(47rm) .T -' 
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and 

(3.9) 


o/_lira ^ 

- / 9 sin{2n6)'y^{—iy J2j{4:Trm) cos{2j6) dO 

“'O 1=1 

y^(— 1 )-^ J 2 j( 47 rTO) / 0 sin( 2 n 0 ) cos( 2 j 0 ) 

1=1 “'o 


/■’1/2 

+ (—l)"'J 2 n( 47 rTO) J 9 sm{2n9) cos{2n9) d9 


- 2n "V _ J2n(47rm) 

i=i 

j^ri 
OO 

= - ^ J2j (471771) 


1=1 


2n 

1 1 


n+j n-j 


J 2 n (471777) 
2n 


= - E 


d2fc-2n(47rm) 


fc = n+l 
k^2n 


E 

fc = — oc 

fc^O 


J2n-2fe(47rm)^ J2n(47rm) 

2n 


__/ ^ J2fc-2n(47rm) -A J 2 «+ 2 fc (471771) ^ J 2 n- 2 fc (47n7l ) ^ J2„(47rm) 

I Zc u / ^ u 


fc = n + l 
fc^ 2 n 


A:=l 


A;=l 


2n 


E ^/2n+2fc(47rm) J2n-2fc(47rm) 

h ^ h ■ 




fe=i 

k^n 


Together with (13.7p and (13. 8p . this implies 

E °° J 2 „+ 2 fc (471771) - J 2 „_ 2 fc (471771) 

k 

k=l 

which, according to (EH), is equal to P 2 n( 47 rm). 

Lastly since f(x,n) is 1-periodic and is of class C^, its Fourier series converges 
to it. This proves the lemma. □ 

The result on ^ Jy( 47 r? 7 i)|^_ 2 ^, arising as Fourier coefficients of a function in¬ 
volving Chebyshev polynomials, inverse trigonometric and some infinite series, is 
now obtained. 

Lemma 3.4. Let g{y,r,x) be defined in (11.121) . For n £ N and 0 < x < 1, the 
identity 


d 




, dv 

m—1 

(- 1 ) 


cos(27r77ia;) 


471 


-{ (f) fe.-. (I) + 


. f ^ — X 
F sm ' 


U2n-1 


1 — X 


. {2 — X 

■ sm ' 


C^2n-l 


2 — a; 
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4 n+l g{m, n, x) + E g{m,n, 1 - x)j 

^ m—1 m—1 ' 


holds. 
Proof. Let 


hiM = (I) C/2„-i (I) +sin-i 


+ sin 


1 — X 


U 2 n -1 


1 — X 


sin 


2 — X 


U2n-1 


2 — X 


(-lln+l / °° °° \ 

h 2 {x, n) = ^ g{m, n, x) + ^ g{m, n, 1 - x) j , 

^ m—1 m—1 ' 

and let 

h{x, n) = hi{x, n) + h 2 {x, n). 

Note that h{l — x,n) = h{x,n). The periodization of h{x,n), as a function of x, 
based on its values in [0,1) makes it an even function of x. Hence its Fourier series 
is given by 

OO 

bo + ^ b( cos(27r£x), 

i=i 

where 

bo = h{x, n) dx, bi = 2 j h{x, n) cos{ 2 ^T^x) dx, 

Jo Jo 

ior £ > 1. The coefficient bo is computed first. Observe that the change of variable 
X = 2 cos 9 yields 

P 1 pTT j 2 

J U 2 n-i (I) sin“^ dx = 2 J sin( 2 n 6 ») d 0 . 


Computing the other three integrals in the similar way, it is seen that 

\ (- 1 )^ 

- / I TT ~ ^ ) sin( 2 nd) dO = 

T Jo '^2 / 


(3.10) 
Note that 


ri (-!)» /-/2 

fii{x, n) dx = 


4n 


rl 2( _/*! ^_^\n+l rl 

h 2 {x,n)dx= ^ g(m,n,x)dx = ^ g{m,n,x) 


dx. 


where the interchange of the order of summation and integration is valid because 
of absolute convergence. 

Let y = m + x—l in the above integral and observe that (m — l + a;)(m + 3 + a:) = 
(m + 1 + x)^ — 4 so as to have 

(_ 1)„+1 ^ p (y + 2 -V(^ + 2)2-4)"" 


pi °° r 

(3.11) / h 2 {x,n)dx= ’ / 

Jo ^ m=l-^™-l 


V(J/ + 2)2-4 
(_l)n+l /■“ (y + 2- v/( 2 / + 2)2-4)'" 


■dy 


22 n+l 


\/( 2 / + 2)2-4 


dy. 
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The change of variable y + 2 = sec 9 transforms the above integral into 

[ + ~ dy = 22"/"^ (sec 0 - tan 6»)2" sec 9 d9 

Jo x/( 2 / + 2)2-4 Jo 


= 2 


2n 


±2n — l 


dt 


Jo 

22n-l 


by another substitution sec0 — tan0 = t. Along with p.lOl) and (I3.11L this yields 
(3.12) ho = 0. 

Next, ht is computed. Using exactly the same approach as in (lT7l) . it is seen that 

2(-i)" 

Jo 


hi{x, n) cos(27r£x) dx = 


TT 


~ sin(2n0) cos (dirt' cos 9) d9 


2 

TT Jo 


0sin(2n0) cos (dir^sind) d9. 
Jo 

Using the formula [101 p. 226] 

OO 

cos(zsind) = Jq(z) + 2^ J 2 k{z) cos(2fc0), 


k=l 


valid for z,9 £ C, and performing a calculation similar to that in p.9l) . one sees 
that 

cl ^ 

(3.13) 2/ hi{x,n) cos(27r£cc) dx = - - —-—(J 2 n+ 2 fe( 47 r^) - J 2 n- 2 fc( 47 r£)). 

“'0 ^ fc=l ^ 

The difficult task now is to evaluate the corresponding integral involving h 2 {x,n). 
This integral is first written in a convenient form using an approach similar to that 
on the previous page, namely, 

2 / h 2 ix,n) cos{2Tr£x) dx 


(3.14) 


^_^^n+1 

/ g(™, n, x) cos(27r£x) dx 

^ m—1 


22r, 


(_l)n+i (y + 2- y'(y + 2)2-4)"”cos(2^£y) 


22ri 


dy, 


\/( 2 / + 2)"-4 

|.7I-/2 

= (_l)"+i / (sec0 — tand)^” sec 0 cos( 47 r£sec( 0 )) c?0 
Jo 

pOO 

= (_l)"+i / cos ( 47 r£cosh(/)) 


where the substitution sec 9 — tan 9 = e was used in the last step. 
C. V. Coates [3l p. 260] showed that the integral 


(3.15) 


pCX) 

(_l)"+i / e“^"''^cos(ucosh(()) dcj) 

Jo 
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satisfies the non-homogeneous second-order linear differential equation 

2n(—1)" cos(m) 


(3.16) 


(Py 1 dy 

du^ M du ^ 


1 - 


An integral whose real part is equal to above integral also appears in Whipple’s 
work [551 P- 106] on diffraction theory. Yet no one seems to have explicitly evaluated 
it. Asatryan [^ Equation (A 9)] has obtained an evaluation of this integral in terms 
of a nested finite sum involving factorials. The following new evaluation of (13.151) 
proved below for m > 0 and n S N, namely 


(3.17) (-1) 


n+1 


e cos (u cosh cf)) dcj) 


= (log (f) - + 1)) J2n{u) 


1 (— 1 )^ 

“ 9 X] - h - {J2n^2k{y^) + J2n-2k{y^)) ” ^ '/2n+2fc ('^^) 






will be useful in completing the proof of this lemma. 

To establish this identity, the existence and uniqueness theorem for second-order 
linear differential equations is employed. Let w{u, n) denote the right-hand side of 
(13.171) . By brute force it can be verified (although it is fairly tedious) that 


(3.18) 


d'^w 1 dw f 4n^ \ 

1 4 °° 

~ ('^2n—1('^) '^2n+l('^)) ^ ^ ^ ( 1) ^J2n-\-2k{'^) 




+ 2n)J 2 n+ 2 fc(xi) + - 2n)J 2 n- 2 fc 



The standard formulas [20l p. 222, 10.6.1] 

and 

(3.19) J,._i(z)-I-J,.+i(z) = — Ji,(z), 

z 

are used in this verification. Substituting j = k — 2n in the sum on the extreme 
right in (]3.18p and simplifying, the right-hand side of (|3.18p is seen to be equal to 


(3.20) — (J2n-l('^) — J2n-\-l{u)) 


^(—1)^(2A: -h 2 n)J2n+2k{u) — 2n ^(—1)^ J 2 n+ 2 fc(?^) 






^ y OO OO N 

~ ^ + 2n) J 2 n+ 2 fc(u) + {-iyjj- 2 j- 2 niu)j 

^ k—1 j—l — 2n ^ 

1 4 

= - (J 2 „-l(u) - J 2 n+iiu)) -J y'(-l)''(2A: -P 2 n)J2n+2k{u) 

k^l 
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d _ . 0 _ 

+ 4 V(-l)" J2n+2fc(u) + - E (-l)"fc./2„-2fe(U- 

/C=l 

Now note that from [El p. 270], 

OO 

(3.21) ^(-1)^(2A: + 2n)J2n+2fc(w) = ^J 2 n-l{u) - 2nJ2n(u), 

k^l 

whereas 

OO 

(3.22) E(~l)^'^2n+2fc('«) = M^2n{u,u) — J2n{u), 

k=l 

where A^{s,z) is the Lommel function of two variables defined by EH p. 537, 
16.5(5)] 0 

OO , n 

—\ / g\y-\-2m 

M,.(S,Z)= E(-ir(-) J.+2m{z). 

771—0 

Further, 

(3.23) 


2n-l 


2n-l 


E i.-^fkJ2n-2k{u) = E - 2fc)J2n-2fe(U 






271-1 


+ ^(—l)^-^2n-2/c('ti) + ^ (—l)^</2n-2fc(w) 




fe = 71 + l 


From [El p. 383, (57.4.15)], 


2n-l 


(3.24) E (“1)^(2’^ - 2k)J2n-2k{u) = - (J2„+i(m) - J-2n+l{u)) - 2nJ2n{u). 


k=l 


Also from [El p. 379, (57.1.25)], 


(3.25) 


fc=i 


E(-l)^'^ 2 n- 2 fe(u) = -^(C 0 SU+ Jo{u)) - M2n{u,u), 


whereas using formula (57.1.19) from [El p. 378], 
(3.26) 

271—1 71 — 1 

E {-^fj2n-2k{u) = (-l)-^(-l)^J2,(u) 

fe = 71 + l j — 1 


= (-l)"+Vo(u) + 


(- 1 )^ 


(cosu + Jo(u)) - M2n{u, u). 


Substituting (13.241) . (I3.25|) and (13.261) in (13.231) yields 


2n-l 


(3.27) E i-^)'^kJ 2 n- 2 k{u) = n(-l)"cOSU + nJ 2 n{u) 


k=l 


- -(J 2 „+l(u) + J 2 n-l{u)) - 2nM2niu,u). 


^The conventional notation Un(s,z) is avoided so as to not get confused with the one for the 
Chebyshev polynomial Un{x)- 
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Substituting (I3.2ip . (13.221) and (13.2711 in (I3.20|) . and using (13.1811 leads to 


--h w 1 - 


<Pw ^ 1 dw 
du^ u du 
2n(—1)" cos u 


4n^ 


QtIj 3 

H- lTJ2n{u) — —{J2n-l{u) + J2n+l{u)) 


2 u 


2n(—1)” cos 1 


where the last step follows from (13.191) . 

Thus the right-hand side of (I3.17P also satisfies the differential equation (13.1611 . 
It is easy to see that both sides of (13.171) agree with each other at u = 0, and their 
derivatives agree at u = 0 as well. The differential equation has a regular singular 
point at a; = 0 with indices ±n. The standard structure of the space of solutions 
shows that if two solutions yi, y 2 have matching values and derivatives at 0, then 
yi = 2 / 2 - This implies that (13.1711 is true for all u > 0 and n G N. 

Now substitute (13.171) with u = diri in (I3.14|) . and then combine the resultant 
with (13.1311 thereby obtaining 


(3.28) be = (log (27r^) - i/;(2n -h 1)) T2n(47r£) - ^ ^ fc(fc+^2n)^^^ 




2k-\-2n 


(dTrt') 




u—2n 


the last step resulting from p.6l) . Since h{x, n) is 1-periodic and is of class its 
Fourier series converges to it. This fact, along with the expressions for &o and he in 
(I3.12|l and (13.2811 complete the proof of Lemma 13.41 □ 


Having proved the above lemmas, we are now ready to complete the proof of 
Theorem 1 1.2 1 Substitute the results of Lemmas 13.31 and 13^ in Lemma and use 
the elementary identity sin“^(x) -|- cos“^(a:) = | to obtain 

(3.29) 

OO 

A{n,x) = (-l)"7r E b2n(47rm) cos(27r?7ix) 

m—1 

(_l)n+l 1 / ~ ^ ^ ^ ^ 

+ 2n ^ 3 {m,n,x) + 2^g(m,n,l-x) 

^ m—1 m—1 

- \ {U2n-l (I) + U2U-1 + U2U-1 + U2U-1 (^) ) ' 

Theorem o now follows from substituting the above representation of A{n, x) in 
(EH) and then making use of the fact that U 2 n-i{—y) = —U 2 n-i{y)- 


4. Recovering Zagier’s formula 

This section shows that Zagier’s exact formula (HH follows by letting a; —>■ 1 
on both sides of the result of Theorem o As is clear from (11.131) . the series 

OO 

(—l)"7r y2Tt(47rTO) diverges, and so the passage to the limit requires some care. 

771=1 
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From (12.Ill , we easily get the big-0 bound 

^(-l)"'7ry2n(47rTO) + cos(27rTOa;) = 0„ . 

This suggests replacing the term (—l)"'7rF2n(47rTO) by (—l)"7ry2n(47rTO) + l/(2-y/m) 
in the statement of Theorem 11.21 The uniform convergence of the series 


^(-l)”7rr2„(47rm) + cos(27rma:) 

m—1 ^ V / 

for X S (0,1) will then allow us to move the limit inside the sum. The consequences 
of adding the term 1/(2-^771) are examined first. Note that 


(4.1) (-1)", y; T2n(47rm) cos(27r77ix) 


OO / ^ \ oo 

= f (-l)”7i’42n{47rm) + j cos(27rma:) - ■ 

m=l ^ V / 


cos{2'!Tmx) 
2^/rn 


An appropriate modification of the second series on the right is obtained from 
the Hurwitz zeta function and its representation [U p. 257, Theorem 12.6] 


C(1-A:.) = ® (cos 


OO 

(f)i: 


cos(27rma:) 


• sin 




oo 

(f)i: 


sin(27 




valid for 0 < a; < 1 and Re(s) > 1. This expansion is also valid for Re(s) > 0 
provided x ^ The special case s = | is used to obtain 

^ cos(27rTOa;) sin(27rma;) 

C(ia:)=> - ^ ^ ^ + > — ^ ^ ■ 

^ ^ ^ y/m ypm 

m—1 m—1 ^ 


Thus 


OO oo / ^ \ 

(—l)"7r y2n(47rm) cos(27rmai) = y^ ( (—l)"7ry2Tt(47rTO) H- -= ) cosi2'Kmx) 

1 1 \ 2 y/mJ 

m—1 m—1 ^ V / 

1 > / I \ 1 sin(27rma:) 

-5<(id+2i: 


The last series is simplified using the identity m Equation (9)] 

OO ^ - OO .- 


m—1 


{m + \/m^ — x'^Y — x^’ 


m—1 


valid for Re(u) > 0 and 0 < a; < 1, and which may be established by applying the 
Poisson summation formula (12.41) to the function J,y{2Trx\t\). Now use the particular 
case u = i and the fact [HI p- 924, entry 8.464.1] that 


to obtain 
(4.2) 


E 


sin(27rmx) 



>^ 1 / 2 ( 2 ) 



X V 

“ 77o 


'^ra=l\l m + \/— x'^ \/to^ — x'^ 
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It follows that 


oo oo . ^ . 

(—l)"7r ^ l2n(47rm) cos(27r77ix) = ^ ( (—l)”7rY2n(47rm) + ^ ^ j cos(27rmai) 

m=l m=l ^ V™/ 


2^^2.^) + 4^ 


E 


1 


m=i \^m + y/m? — x'^ \/m? — x"^ 

The result of Theorem fOl is now expressed as 

^ / 1 \ 11 

B^nix) = ^ f (-l)"7ry2n(47rTO) + j cos(27rTOa;) --C(i,a;) + ^ 

m—\ ^ V / V 

^ X 1 (3 — a; — Y^(l — x)(5 — x))^" 

[2x/2v'l + \/l - xWl - 22 "+i ^(l-a;)(5-a:) 

1 


E 


2^/2 m=i Y^TO + 1 + y^(m + 1)2 — a;^' -\/(m + 1)^ — x^ 

+ l (U2n-1 + U2n-1 (|) + U2n-l + U2n-1 (^)) 

^ / OO oo \ 

~^ 92 n+l E g{m, n, x) + ^ g{m + 1, n, 1 - x) . 


\m—l 


m—1 


Now let X t 1, use the facts that C/2rt-i(0) = 0, t/ 2 n-i(l) = 2n and that U 2 n-i{x) 
is an odd function of x, along with the value 


lim 


1 (3 — X — •\/(T^^'x)(5^^^)2 


1- V2V2Vl + Vl-x2^/r^ 22 "+i v'(1-x)(5-x) 

1 


2 4\/2’ 

to obtain 


BIM) = f:((-ir’'L-«™') + 5^)-lc(5) + j 

m—1 ^ V / \ / 


y, ^ 

m=i y/m{m + 2) W m + 1 + y^m{m + 2) 


4 ^ 2^2 


+ E 


^Jm{m + 4) 


Vm + 4 — ■v/m^ 

2 j 


where the identity m + 2 — ^m{m + 4) = ^{y/m + 4 — i/ro)^ was used in the 
simplification of the resulting series on the extreme right side. The last step is to 
use 

2^ _ 2 

m + 1 — m{m + 2) = - (-v/m + 2 — v^) 
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to evaluate, as a telescoping series, 


m=i y/m{m + 2) Wto + 1 + m{m + 2) 


^+1 


Finally use the result + R [3 Formula (10.16)], [SOj p. 5] to obtain 

Zagier’s formula (ini). 


5. Periodicity and the odd-index case 

The original motivation that led us to the study of the modified Bernoulli num¬ 
bers was the curious phenomenon that : n G N} is a periodic sequence. 

Since the proof of Theorem 11.31 is similar to that of Theorem 11.21 it is omitted. 
However, a new derivation of the periodicity of {B^^+i} using Theorem ll.3l is given 
below. 


Corollary 5.1. For n G N, 

(- 1 )” , 1 


B. 


2n+l 


2 V 2 


4 V 3 


Proof. The series in Theorem 11.31 may be written as 

OO 

(-ir-E y2n+i(47rm) sin(27rTOa;) 

m—1 

°° / 1 \ 1 °° 
= E f (-l)” 7 '' 42 n+i( 47 rm) -b j sin(27rma:) - 2 E 

m—1 ^ V / m—1 

/ 1 \ 1 

= E (-l)"7r42n+i(47rm) -b sin(27rma:) - 

00 

E 


sin(27rTOa;) 


m=i \/m F \/rrF- — 'Jvnf — 

using (14.21) . Now replace the above representation in the formula from Theorem ll.3l 
to obtain 


B. 


2n+l 


(x) = E f(-l)"7rF2n+i(47rm) -b sin(27rTOx) 

771=1 ^ V / 


E 


2’\/2^^i \/m F \/rn^ — \/rvF — 


^1 ^ (2 -b X — y^x(4 -b x))^"+^ 1 

4 y 22"y/x(4-bx) 


X — 1 


C/2„ 


X — 2 


-b 


1 


22n+2 


E 5 (w -b 1 , n -b i, x) - E 5 (w, n -b 5,1 - x) 


V 771=1 


771=1 
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To obtain the final result, let x —>■ 0 in the identity above. The use of the dominated 
convergence theorem, the evaluation 


lim 


1 {2 + X — yj x(4 + 1 


0 \ 22'’ 


= —2n — 1, 


y^x(4 + x) \fx ^ 

the parity of the Chebyshev polynomials and the special values 


b' 2 n( 0 ) = (-1)”, 172n(l) = 2n + 1, C/ 2 „ (i) = sin 


{2n + 1)71 


complete the proof of the stated formula. 

6. Asymptotics of Zagier polynomials 


□ 


Corollarv ll.dl is proved here. The proof is given only for the even-index case, the 
one for odd being similar. 

First assume x ^ The idea is straightforward: we divide both sides of 

Theorem 11.21 bv (—l)"7ry2n(47r) cos(27rx) and show that the resulting right-hand 
side approaches 1 as n —>■ oo. However, in the case of the Bessel function series, 
justification of the interchange of the order of limit and summation is needed, and 
which is interesting in its own right. To that end, it is first shown that 


( 6 . 1 ) 


lim 


1 


CX) 


— 1)” 711271 (47rm) cos (27rmx) = 0. 


rt^oo (—1)” 7rl2n (47r) COS (27rx) ^ 

Indeed, use dSI) to rewrite the above limit in the form 

(^(-1)” 7ry2r7 (47rm) -f cos (27rmx) 


lim y 

n—^oo 

m—2 


(—1) TTY2n (47r) COS (27rx) 


— lim 


-.->■00 (—1)" 7rl27i (47r) COS (27rx) 

The latter limit is equal to zero since (12.21) implies that 


cos [2'Kmx) 
2y/m 


( 6 . 2 ) 


lim 


= 0 , 


n->oo y2n(47r) 

and also because the series X]m =2 converges for 0 < x < 1. 

It is now shown that the order of limit and summation can be interchanged in the 
former limit. This requires the hypotheses of Lebesgue’s dominated convergence 
theorem for series to hold, namely that 


(6.3) 


lim 

n—foo 


((-!)” ^^277 (47rm)+2^) 


= 0 , 


(-l)”7rr2„ (47r) 
and that the sequence {s„(to, x)}))T]^ , where 

(^(-1)” 7ry277 (47rm) -f 2 ;^) cos (27rmx) 


(6.4) 


Sn(m,x) = 


(—1)" 71^277 (47r) cos (27rx) 
is uniformly bounded in n. To prove (16.31) . note that (12.21) implies that 

1^277 (47rm) 
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as m > 2. Together with (I6.2L this proves (I6.3p . 

In order to prove (I6.4L we first use the fact [9l p. 927], [27l p. 446] that for a; > 0, 
the function 

(x) + (a:)] 

decreases monotonically, if u > i. In particular, for m,n gN, 


iirm [ J|„ (dirm) + (dirm)] < dir [ j|„ (dir) + Y^^ (dir)] . 

From p. 227, formula 10.14.1], we have j (x) j < 1 for u > 0, a; £ R, so that 

(6-5) Y^^ (dirm) < j|„ (dirm) + Y^^ (dirm) < — [ J|„ (dir) + Y^^ (dir)] 

< — [1 + (dir)] . 

m 

Next, it is shown that for n £ N, 


for some absolute constant c. Note that the asymptotic formula (12.21) implies that 
the sequence {y2n^(47r)}^_^ tends to 0 as n —>• oo. So for n sufficiently large, say 
n> no, the inequality y,; 1 . . < 1 holds. 

If we can further show that T2„(47r) ^ 0 for any n < no, then (16.61) will be 
proved as one can take c to be max (l, F 2 (dir), 14 ( 471 ), • • • ,y'2(no-i)(4'T)). To show 
that this is indeed true is the objective of the following lemma. We could not find 
a reference to it in the literature, and our proof of it is short and nice, hence given 
here. 


Lemma 6.1. For any positive integer n, T 2 n( 47 r) ^ 0. 

Proof. Let denote the fc**' zero of the Bessel function W {x). The last line on 
page 68 in [ 8 ] implies the inequality 

> u + fcTT - i (i/>4,fc£N). 

Thus the first zero of F 2 „(x) satisfies y 2 n,i > dir if 2n + tt — 4 > dir which happens 
when n > ^ + I Ri 4.96239. Hence F 2 n (dir) 7 ^ 0 for n > 5. Also, it can be checked 
that r 2 n (47r) 74 0 for 1 < n < 4, since Y 2 (dir) r; 0.134559, 1^4 (dir) Ri -0.0357975, 
Ig (dir) Ri —0.14694 and Is ( 47 r) Ri 0.246447. This completes the proof of the 
lemma. □ 


The result in Lemma 16.11 along with the previous discussion now proves (16.61) . 
Now divide both sides of (|6.5I) by T2^(47r) and use (16.6|) to obtain 

l2n(47rTO) ^ \/c + 1 


m 


l2„(47r) 

thereby proving dSH). 

Thus, from Lebesgue’s dominated convergence theorem, it is seen that 


lim y 

T3.—^ 


m—2 


(^(-1)" 7rl2n (dirm) + cos (27rma:) 

(—1)” 7rF2n (dir) cos (27ra;) 
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E lim 

n—>00 


(^(-1)" TTY 2 n ( 47 rm) + 2 ^) cos {2Trmx) 
(— 1 )" 7 rF 2 n ( 47 r) cos {2ttx) 


= 0 , 


as can be seen from (lO) and (lOl) . This proves (EH). 

Our next task is to show that 
(6.7) 

g (m, n, 1 — x)^ = 0 . 

Again using the series version of Lebesgue’s dominated convergence theorem, we 
observe that as n —>■ 00 , 


1iTn - 

n^oo 22"+1 (—I)" 7rl2ri (47r) COS (27rx) 


f 00 00 

^ g{m,n,x) + ^ 

\m—l m—l 


^ g (to, n, x) - 5 (1, n, x) = 


^4n 


^2 + X + ^x (x + 4)^ -^x (x + 4) 


so that by ( 12 . 21 ) . 


lim 


Y.Z=i9ijn,n,x) 


rt-J-oo 22"+1 ( — 1) 711271 (47r) COS (27rx) 


1 (-l)"+^22"-iy7^ 

- — lim - 

cos (2tx) x/^a 7 (x + 4) (2 + x + y^x (x + 4 )) 


•^71 ^e7r^2n 


2 n 


= 0 . 


Similarly, 


1 


lim - 

n-700 22 "+i (—1)” 711271 (47r) COS (27rx) 


Y.9^'' 

m—l 


, 1 — x) = 0. 


Thus (16.7p is proved. Finally the fact that 
( 6 . 8 ) 


(^277-1 (^) + f/277-1 (f) + ^277-1 (^) + f/277-1 (^)) 

U.-700 4 (—1)" 7rl277 (47r) cos (27rx) 


easily follows from (|2.3I) and (|2.2I) . Thus (16.11) . (16.71) and (16.81) along with Theorem 
Oprove (frra for x ^ i ^ 


4; 4- 


When X = j,j, the first term of the Bessel function series in Theorem O is 
zero. However, a logic exactly similar to the one above can be worked out starting 
with the second term of the series, namely, (—l)”“''^ 7 r 1277 ( 871 ). 


7. Proof of the Zagier-type formula 

This section is devoted to proving Theorem 11.51 Let x = 1/2 in Theorem 11.21 
and note that 1 / 277-1 (x) is an odd function of x. This gives 

(7.1) 


BL 


= ^(-l)-+"Ty277(47rTO) + ^^ 


m—l 


(to + I - 

("* + !) 
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Now let X = —3/2, k = 2, and replace n by 2n in (|1.14p . Along with ()7.ip . this 
gives 

1 (V 2 - 1 ) , n 


(- 1 ) = 

^ / m—1 ^ V / 


- 2 1 U2n-1 


+ U2n-1 


(^ +1 - - i) (^ + i)) 

^ m=l \l {^~ \) I) 

where we also used the well-known identity 


m—1 




valid for Re(s) > 0. 

Add the corresponding sides of the above equation to those of Zagier’s formula 
(ini) to obtain 

BL = E (1 + (-I)’”) (^(-l)"^r2„(47rm) + 


1 


— n — — \ U2n-1 T + U2n-1 




—i y 

22n I / ^ 


+ 2— ^Jm{m + 4)^ 
y m{m + 4) 


2n 


(7.2) 

Note that 


^ (w + I - v^(m - ^) (m + 1)/ 

^ V-l) (-+1) 


2n 


E (l+(-l)™) ((-l)”^r2„(47rm) + = 2 E ((-l)"^F2„(87rm) + , 

m—1 ^ V / m=l ^ V / 

and that the sum of the two series in (IL21) can be written as 

CXD 

E 

m—1 

^ °° + 4 — ^m{m + 8)^ 

“ 92 n-l 


f (¥ + 2-/2fi(2-+4)) ” 

(2^+2-/ 

V 2n 

(2^4) (2^4+ 4)) 1 

1 /2-(2^+4) 


1)(2™-1+4) J 


™=i y'm{m + 8) 

This completes the proof. 


8. Some open problems 


We conclude this paper by discussing three open problems. 
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8.1. Proving Lemmas 13.31 and 13.41 throuerh Poisson summation formula. 

We begin with a lemma. 

Lemma 8.1. For Re(i 2 ) > 0 and 0 < a; < 1, the identity 

OO 


m—1 


J^{iTrm) cos{2TTmx) 


cos (z^arcsin (l)) cos (z^ arcsin (^^)) 

(47r)2 — (27ra;)2 (iTr)^ — (27r(a:: + 1))^ 

^ cos (z^arcsin (h^)) ^ cos (z/arcsin (^^)) 


y^(47r)2 - (2 t( 1 - a;))2 y^(47r)2 - (2 t( 2 - a;))2 

(4^)-sin (if) 


-E 


-E 


{m + x))2 — (47r)2 ^27r( m -\- X ) +m + x)Y — (47r)2^ 
(471)'^ sin (if) 

(m — x))2 — (47r)2 ^27r(m — a:) + y^(27r( m — a:))2 — (47r)2^ 


holds. 


Proof. Apply the Poisson summation formula (12.41) with f{t) = J,y(47r|t|) cos(27ra;|t|), 
Re Z7 > 0, and use the integral evaluation Up. 717, 6.671.2] 0 



Ju{at) cos(/3t) dt 


cos( 1 / arcsin 




(0 < /3 < o), 


, , - \ i XI. , (0 < a < /3). 


□ 


8.1.1. An alternative approach to proving Lemma \3.3[ Consider the double series 

OO OO ^ 

(8.1) 3(a:, ■«) := E E ('^2n+2fc(47rTO) - J2n-2fe(47rm)) cos(27rma;). 

k—1 m—1 

We wish to use Lemma [Q twice to simplify the result. Observe that it requires 
Re(z/) > 0, so this requires to treat the case k = n separately. In this situation, 
entry 8.522.1 in [9] states that 

OO 

(8-2) E Jo{mu) cos{muv) 

m—1 

_ 1 1 1 1 

2 \/y^ — (Stt^ + uv)^ uVl — v"^ ^ \/'uf — (27rf — uv)^ ’ 

where rz > 0, 0 < u < 1, 27rj < u{l — v) < 2{j + l)7r, 2d'iT < u{l + v) < 2{d + l)7r 
and j + 1, d + 1 G N. Now use Lemma 18.11 twice in m, cLiid. clIso 11 S 6 (|8.2jl to 


■^The condition that a and /3 be positive is missing in the reference. But it can be found, for 
example, in [27l p. 398, Section 13.4]. 
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obtain after some simplification 

, 1 1 “ 1 rsin(2nsin-1 (f)) sin (2fcsin-1 (I)) 


fc=l 


\/4 — 


The definition in 


now yield 
3{x,n) = ^ 


sin(2nsin i (^^)) sin (2fcsin ^ 

— (a; + 1)^ 

sin (2nsin-i (^^)) sin (2A;sin-i (^^)) 
V'd - (1 - a;)2 

sin(2nsin-i (2^))sin(2fcsin-i (^^)) 
^ v'4-(2-x)2 

and the classical formula [H p. 46, formula 1.441.1] 
^ sinj0 TT — 9 
1=1 -1 


for 0 < 0 < 2?! 


2n 2Tr 
cos-1 


cos 1 c/2n-l (I) + cos" 
U2n-1 


1 — X 


1 — X 


-1 f x + l 

V 2 
2 — X 


U2n-1 


COS 


U2n- 


X + 1 
2 

2 — X 


Comparing with the result of Lemma 13.31 and the identity (13.41) , it is easily seen 
that Lemma 13.31 can be proved this way if the following problem, which we leave 
for the interested reader, can be solved. 

Problem 1. Prove that for 0 < x < 1 and n £ N, 


OO OO 


1 


^ ^ - (J2„+2fe(47rm) - J2„-2fc(47rTO))cos(27rTOx) 
k—1 m—1 

OD OO ^ 

= ^ ^ - (J2n+2fc(47rTO) - J 2 n- 2 fe( 47 rm))cos( 27 rmx). 


m=l k=l 


8.1.2. An alternative approach to proving Lemma \3.4\ Using Lemma l8.11 it can be 
seen, after some simplification, that 


d 


lim Ji/(47rm) cos(27rTOx) 


v^ln dv 

{-ly 


An 


U2n-1 


sm 


1 — X 


• -1 f x + 1 

-n I 
1 — X 


U2n-1 


X + 1 


sm 


2 — X 


U2n-1 


2 — X 


(- 1 ) 


n+1 


4^+1 


^ g{m, n, x) + ^ g{m, n,l-x) 


If we compare this result with that in Lemma 13.41 it is easily seen that one can 
prove the lemma in this alternative way provided the following interchange of the 
order of summation can be proved. 
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Problem 2. Prove that for Re(i^) > 0, 0 < a: < 1 and n G N, 


°° 5 d 

E —Ji/(47rm) cos(27rTOa;) = lim — > J„{ATTm) cos{2TTmx). 

hij „ V^2n Hi/ ^' 

y—2n 


m—l ' L/—zn m—1 

This, too, is left as an open problem for the reader to prove. 

Remark. Problems analogous to above can be formulated for the odd-indexed 
case. 


8.2. A relation bet-ween R* (— §) and R*? Theorem 11.51 arose while trying to 
find a relation between R* (—|) and R* similar to the relation [22l p. 4] 

Bn Q) = (2'-” - 1) R„ 

that exists for Bernoulli polynomials. Note that while the result Bn{l — x) = 
(—l)"R„(x) shows symmetry along x = 1/2, it was established in [71 Theorem 
11.1] that 

(8.3) B:{-x-3) = i-irB:ix), 

so that the symmetry for Zagier polynomials is along x = —3/2. Since the theory of 
Zagier polynomials nicely parallels that of the Bernoulli polynomials, it is reasonable 
to look for a relation between R* (—|) and R*. Of course, it is clear from (E3D 
that R 2 „+i (~|) = 0. On the other hand, the sequence {RJn+i} is 6-periodic and 
takes the values -j, j, j, -j}, readily implying 

^2n+l (-0 = • ^2n+l- 

Hence only the case of even indices is of real interest. We have not been able to 
find such a relation yet, if at all it exists. 
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